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Abstract 



The interest in the study of spin-Peierls phenomenon has increased 
with the discovery of copper germanium oxide (CuGeO^) as an example of 
spin-Peierls systems. The availabihty of this inorganic compound as large, high 
quality single crystals, allow experimental measurements which aid in under- 
standing the spin-Peierls phenomenon. This report details the numerical and 
analytical work done on the various magnetic phases that occur in a quasi one- 
dimensional spin-Peierls system (like that of copper germanium oxide). Taking 
the system to be finite, the bond displacements, single particle energies, energy 
eigenfunctions, spin densities and phonon frequencies are calculated numeri- 
cally for all possible magnetizations, which includes both the commensurate 
and incommensurate phases, at zero temperature. Solitons of fractional spin 
are found to be present in the system for appropriate magnetizations. The uni- 
form to dimerized phase transition occurring in this system is also modelled 
numerically. The single particle and ground-state energies arc later obtained 
theoretically for the uniform and dimerized phases. In addition, the lattice spin 
displacements arc analytically modelled for a particular magnetization of the 
system. 
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Chapter 1 



Introduction 



1.1 Spin-Peierls systems 

A.n interest in the theoretical study of quantum spin systems exhibiting the 
spin-Peierls transition has grown ever since the discovery of the first inorganic 
spin-Peierls compound, copper germanium oxide (CuGeO^). A spin-Peierls 
transition is a structural transition associated with quantum spin systems where 
spin phonon coupling leads to distortions of the lattice. This magneto-elastic 
phenomenon gives rise to a variety of phases in the system as a function of 
field H and temperature T This report is aimed at analyzing these various 
phases of the spin-Peierls system and the associated physical quantities. 

In this chapter, the spin-Peierls systems are introduced by first dis- 
cussing the structure of the spin-Peierls compound CuGeO^. Its phase diagram 
is studied later and the section is concluded by writing the model Hamiltonian 
for this spin-Peierls system. 

1.2 Structure of CuGeO^ 

The crystal structure of CuGeOs is shown in Fig. ( |1 . 1|) . The crystal is con- 
structed of Cu02 ribbons connected via GeO^ tetrahedra. The Cv^^ ions in 
3(i^ configuration form antiferromagnetic spin 1/2 chain in c direction. There 
are two different types of shown in the figure in black (02) and grey (01). 
The Ge^~^ ions connect adjacent Gu02 ribbons. The lattice parameters of the 
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Figure 1.1: Structure of the spin-Peierls compound CuGeO^ 




Figure 1.2: Projection of Cu02 ribbons in (a) three dimensions, and (b) a-b 
plane. 

orthorhombic unit cell are a = 4. 81 A, b = 8.43A and c = 2.95 A ||l|. The Cu- 
02-Cu (with a bonding angle t] ~ 99") is responsible for the super exchange 
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intrachain couplings. The lattice distortions that occur in the dimerized phase 
obtained by the neutron diffraction measurements are shown in Fig. ( |1 . 21 ) . The 
arrows indicate the direction of atomic displacements. 



1.3 (H,T) Phase diagram of CuGeO^ 




Figure 1.3: (H,T) phase diagram of the spin-Peierls compound CuGeO^ 

The (H,T) phase diagram of a spin-Peierls compound CuGeO^ which 
closely resembles the theoretical predictions is shown in Fig. ( |1.3| ) The 
diagram consists of three phases: the uniform U phase, the dimerized D phase 
and the incommensurate / phase. At high temperatures, the system is in the 
uniform phase. In this phase, the magnetic chains have uniform spacing / and 
one exchange coupling J between the neighbouring spins. As the temperature 
is decreased, a second order structural transition is observed at Tsp ~ 14k in the 
case of CuGeOs- In zero or low fields, this transition occurs and it corresponds 
to a dimerization of lattice. In that dimerized phase, two lattice parameters li 
and I2 and two exchange couplings Ji and J2 are needed to define the magnetic 
chains. Increasing the magnetic field from the D phase yields a new transition 
occurring at specific value of the field He, ~ 12. 5T at zero temperature in 
CuGeO^. This field induced transition corresponds to a new deformation of 
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the lattice. Above He, the lattice becomes incommensurate i.e., the period of 
modulation of spin density is generally not a rational multiple of the lattice 
constant / as against the uniform and dimerized phases. Eventually, as the 
applied magnetic field becomes very large, all the spins get polarized and we 
get a ferromagnetic state. 

1.4 Modeling the spin-Peierls compound CuGeOs 

The Hamiltonian Hgp used to model the spin-Peierls system is given by 

N 

Hsp = J^[^- K ^n+l - Un )] 
n=l 

N u N 

-hY.si + ^Y: K+1 - u^f + E ^ (1-1) 

n=l ^ n=l n=l ^^^^ 

where is the number of sites, h is the strength of the magnetic field, Sn 
represents the n*^ spin, Un is the displacement of lattice point n, Pn represents 
the momentum of lattice point n, M is the mass of each lattice point and A is 
the positive constant which is a measure of the coupling that exists between the 
lattice points and spins. In the above equation, periodic boundary conditions 
are assumed i.e., Sn+i = Si. 

The first term in the Hamiltonian gives an antiferromagnetic inter- 
action between the nearest neighbour Hciscnberg spins along the c direction of 
the lattice with their exchange coupling J modified to take into account the 
spin-phonon interactions to the first order in bond displacements. Thus J is re- 
placed by J [1 — A( Un+i — Un )]. The interactions between the spins in the other 
two lattice directions and between the spins other than the nearest neighbour in 
c direction of the lattice are neglected. Considering the distance of separation 
between the spins, this is a good approximation and the system can be ana- 
lyzed as a quasi one-dimensional (ID) system. The second term is the Zeeman 
term giving the effect of magnetic field on the system with net magnetization 
Y,n=i ^n- third term is the harmonic energy of the lattice and the final 
term represents the kinetic energy of the lattice. 

This Hamiltonian is simplified later to calculate the ground state 
energies and other physical quantities of interest discussed in the next chapter. 



Chapter 2 



Numerical Results 



2.1 Simplification of the Hamiltonian 

The full Hamiltonian for the general spin-Peierls system is shown in 
the previous chapter, Eq. (|1.1|) . It is difficult to work with this Hamiltonian to 
calculate quantities of interest like single-particle energies and wave functions 
at zero temperature, spin densities etc. Hence the full Hamiltonian is simplified 
using the physical arguments discussed below so that the quantities calculated 
from the simplified Hamiltonian reproduce the physics of the full Hamiltonian 
to a very good approximation. 

2.1.1 Approximations involved and their justification: 

• Soft phonon limit: The time scales involved in the motion of lattice are 
known to be much larger than the time scales involved in the motion 
of spins arising due to the interactions between the spins and phonons. 
Hence the lattice points are initially taken to be at fixed positions while 
calculating the unperturbed single-particle energies at zero temperature 
and zero field. At the end, the lattice movements are however taken into 
account using the perturbation theory and the perturbed single-particle 
energies due to motion of lattice points are calculated to second order in 
bond displacements. 

• XY spins: The generalised Heisenberg spins in the model are replaced 
by the XY spins, i.e., the spin-spin interactions in the z direction are 
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neglected. This is done to simplify the calculations. When the spins 
are transformed into Jordan- Wigner fermions, the interactions in the z 
direction involves four spinless fermion operators [Q. It is known from 
numerical analysis that reducing the Heisenberg spins to XY spins does not 
qualitatively affect the phase diagram of a quasi ID spin-Peierls system. 

The simplified Hamiltonian, Hg, thus obtained using the above ap- 
proximations is given by 



N 
n=l 



N 



N 



+ ^n^n+l + Tj'Yl ^n+^ ("^-^^ 

n=l n=l 



where ^n^\ — ~^ (""n+i — Wn) is the bond displacement and A is the positive 
constant which is a measure of the coupling that exists between the lattice 
points and spins. 

It is to be noted that in addition to the mentioned approximations, 
a unitary transformation of the form shown in Eqs. ( |2.2| ) is performed on the 
Hamiltonian given in Eq. ( |1 . 1| ) to reverse the sign of the first term. 



/ 1 "i" QX 

SI ^ {-lysi 

~^ (2-2) 



2.1.2 Jordan- Wigner transformation 

Jordan- Wigner transformation transforms a spin Hamiltonian into a spinless 
fermion Hamiltonian p|, ^, Here, two spinless fermion operators, the cre- 
ation and annihilation operators, are defined as shown below 



ct 



n-1 
n-1 

n 



(2.3) 
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where 

c- _ qx _ ■ qy 

q+ _ qx \ o qy 

and (j's are the Pauh matrices. Using the properties of Pauh matrices, it follows 
from the relations (|2.3| ) between spin operators and fermion operators that the 
total magnetization, m, of the system is 



n=N 

m=Y.S:, = Nj.-- (2.4) 

n=l ^ 



where Np = J2n=i C^n^n is the number of Jordan- Wigner fermions. 

Using the Jordan- Wigner transformation on the simplified spin Hamil- 
tonian ( ^.11 ), a spinless fermion Hamiltonian is obtained as shown in Eq. 



H 



J 
2 



.n=l 



N 



n=l 



N 



^ n=l 



(2.5) 



2.2 Single particle energies 

The Hamiltonian used to calculate the single-particle energies at zero tempera- 
tures for all phases ( for all possible magnetizations) is given by 



E 



n=l 
N 



2 



71+7 



CnCn+1 + C^+lCn 



n=l 



n=l 



N 



C^Ci + C^Cn 
(2.6) 



It can be seen that the boundary term of this Hamiltonian has a 
relative minus sign compared to the Hamiltonian shown in the Eq. ( p.5| ). This 
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appears to break the translation invariance. However, Eq. ( p .61) can be made 
translation invariant by transforming all the fermion operators as 



Cl ^ e--/^C+ (2.7) 

which distributes the minus sign at the boundary over all the sites equally. As 
— s> oo, such a transformation produces a vanishingly small change in the 
single-particle energies. 

The system size N is chosen to be 100 spins. The constant c in the 



final term in the Hamiltonian ( p.6|) showing the harmonic energy of the lattice 
is set equal to unity. The lattice constant / is also fixed to be unity so that 
all lengths which are multiples of the lattice constant are now just integers. In 
order to obtain the single-particle energies for any particular magnetization or 
any particular filling of Jordan- Wigner fermions, Np^ of the ground-state, it is 
necessary to obtain the set of bond displacements, 5„+i's which minimize the 
Hamiltonian given in Eq. (|2.6|) . The procedure followed to achieve this is 



Set all (5„+i's to zero initially. 



• Increment or decrement each ^n+\ separately in steps of 0.0005 and change 
each of the remaining spins by 0.0005/99 with the reverse sign, so that for 
each such change the sum of all displacements continues to remain zero 
as required for a periodic lattice. 

• Calculate the new ground-state energy for each such change and iterate in 
a loop running over all the lattice sites, so as to minimize the Hamiltonian. 

Then the eigenvalues of the Hamiltonian with these bond displacements, 's 
give the single-particle energies for the particular filling at zero temperature. 
The eigenfunctions similarly give the single-particle wave functions. The bond 
displacements and single-particle energies are obtained for all possible magne- 
tizations, few interesting ones of which are shown in the following sections. 
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2.2.1 For zero magnetization (for Np = 50) 

The bond displacements for the case of zero magnetization, i.e., for Nf=^^-, 
obtained by minimizing the Hamiltonian in Eq. ( |2.6| ) using the above procedure 
are shown in the lower part of Fig. (|2.1| ). It can be seen that the value of the 
displacement of the bond alternates in sign between any two consecutive lattice 
spins with its magnitude remaining constant. This corresponds to dimerization 
of the lattice. The magnitude of each of the bond displacements, 5, is 0.0642. 



A 

I 
I 
I 




_g -I I \ \ \ I \ \ I \ I I 

■ 10 20 30 40 50 60 70 80 90 100 

n > 



Figure 2.1: The bond displacements and expectation values of for zero mag- 
netization. 

The single-particle energy spectrum obtained by diagonalizing this 
Hamiltonian with these bond displacements is shown in Fig. ( p.2|) . It is observed 
that 

• There are an equal number of positive and negative energy levels. 

• For every positive energy level there is a corresponding negative energy 
level which is equal in magnitude but opposite in sign. 
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1.5 




-1.25- 

_.| 5 1 I I \ \ I \ \ I I 

■ 10 20 30 40 50 60 70 80 90 100 

Site number (n) > 



Figure 2.2: The single-particle energy spectrum of the zero magnetization. 

• Each level is doubly degenerate. Exceptions to this fact are the maximum 
level (1), the minimum level (-1), the largest negative level (-0.0642) and 
the smallest positive level (0.0642) which are non-degenerate. 

• There is a gap between the largest negative level and smallest positive 
level which is twice of 5. It is to be noted that all the quoted energies are 
in units of J. 



The ground-state energy of this system is obtained by summing the 
lowest fifty levels i.e., all the states lying below the gap (negative energy states). 
This is the state of the system at zero temperature and its energy is -31.8634. 

The upper half of Fig. ( p.l|) gives the variation of spin density with 
the site index. It can be seen that this remains constant at half throughout the 
chain as expected. Also the bond displacements, single- particle energies and 
spin densities remain independent of an applied magnetic field for small values 
of the field because the net magnetization in this case is zero. 
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2.2.2 For magnetization=l (for Np=h\) 



0.6 r 
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Figure 2.3: The bond displacements and expectation values of for magneti- 
zation=l. 

The bond displacements for the case of 51 fermions i.e., one away 
from dimerization, obtained by minimizing the Hamiltonian in Eq. (|2.6| ) for 
the case of zero magnetic field using the same procedure, are shown in the 



lower part of Fig. (|2.3|) . It can be seen that the displacements of the bonds 
oscillate around the value zero and have an envelope whose periodicity is 50, 
their maximum magnitude being 0.058. The pattern obtained is similar to that 
of beats in sound waves. 

The single-particle energy spectrum obtained in this case is shown 
in Fig. ( p.4|) . The features of the spectrum that get modified compared to 
Np = 50 are 



There is an increase in the number of non-degenerate levels by four. The 
doubly degenerate levels lying next to the non-degenerate levels of the last 
case become non-degenerate. 
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100 



Figure 2.4: The single-particle energy spectrum for magnetization=l. 



• The size of the energy gap existing between the greatest negative level 
and smallest positive level falls from 0.0128 to 0.02. 



The ground-state energy of this system is now obtained by summing 
the lowest fifty one levels and its value is -31.8219. 

The upper half of the Fig. ( p.3|) gives the variation of spin density 
with the site index. The pattern obtained is seen to be similar to that of bond 
displacements except that the spin density maximum occurs at the site where 
the bond displacement is minimum and the oscillations within each envelope 
remain always above the value half. Also, the area under each of these envelopes, 
say from n = 35 to n = 85, is half and the sum of the areas under the two 
envelopes is equal to unity. This implies the extra (51'^*) fermion distributes 
itself into these two envelopes with equal probability. Thus, the solitons with 
half integral spin are seen in the system. 
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2.2.3 For magnetization=49 (for ^¥^=99) 




Figure 2.5: The bond displacements and expectation values of for magneti- 
zation=49. 

The bond displacements for the magnetization of the system being 
49 or fermion filling being just one less than the maximum allowed, for the case 
of zero magnetic field are shown in the lower part of Fig. ( p.5|) . It can be seen 
that the system has a single peak (~ 0.18) in bond displacements localized over 
a small region chosen randomly by the system by spontaneously breaking the 
translation symmetry. 

The upper half of Fig. ( p.3| ) gives the variation of spin density with 
the site index. It can be seen that there is a dip in the spin density over the 
region where there is a peak in bond displacements, showing that there is a 
localization of the missing (lOO*'^) fermion over a small region. 
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"'O 10 20 30 40 50 60 70 80 90 100 

n > 

Figure 2.6: The bond displacements and expectation values of for magneti- 
zation=18. 

2.2.4 For N=102, magnetization=18 (for iV^ =69) 

The bond displacements for a system size of 102 spins, Np being one away 
from trimerization i.e, 69, obtained by minimizing the Hamiltonian in Eq. ( |2.6| ) 
for the case of zero magnetic field are shown in the lower part of Fig. ( |2.6| ). 
The bond displacements of the lattice spins are seen to oscillate around the 
value zero, their maximum magnitude being 0.022 and have an envelope whose 
periodicity is 34. The pattern is similar to Np = 51 case. Also, the single- 
particle energy spectrum is seen to be gapless in this case. 

The upper half of Fig. ( p.3|) gives the variation of spin density with 
the site index. The pattern obtained is similar to that of bond displacements ex- 
cept that the spin density maximum occurs at the site where the bond displace- 
ment is minimum and oscillations within each envelope remain always above 
the value 0.667. Also, the area under each of these envelopes, say from n = 24 
to n = 57, is 0.333 and sum of the areas under all three envelopes is unity. This 
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means that the extra ( 69*'') fermion distributes itself into the three envelopes 
with equal probability. Thus, solitons with fractional (1/3) spin are seen here 
in the system. 



2.3 Variation of ground-state magnetization with 
magnetic field 




_35 1 \ \ \ \ I I I I I I 

-50 -40 -30 -20 -10 10 20 30 40 50 

magnetisation (m) — > 

Figure 2.7: The variation of ground-state energy with magnetization m for zero 
magnetic field h. 

The single-particle energy spectrum, the ground-state energies and 
the spin densities at zero temperatures are calculated for each of the magne- 
tizations (0 to 50) possible in the system. Fig. ( PTT] ) shows the variation of 
ground-state energy Eq with magnetization m, with the magnetic field term set 
to zero in the Hamiltonian (|2.6|) . 



The variation of ground-state magnetization of the spin-Peierls sys- 
tem with magnetic field at zero temperature is now studied numerically, by 
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-1.5 -1 -0.5 0.5 1 1.5 

Magnetic field (h) — > 

Figure 2.8: The variation of ground-state magnetization m with magnetic 
field h. 

varying the magnetic field h from -1.1 to 1.1 in steps of 0.001 in the Zeeman 
term, h Z]^=i[Cn C*n — |], of the Hamiltonian ( p.6|) . The ground-state energies for 
each of the possible magnetizations are calculated and the magnetization which 
minimizes the ground-state energy is identified for each value of the magnetic 
field. The plot in Fig. (|2.8|) shows the variation of ground-state magnetization 
m with magnetic field h. It is seen that the saturation magnetic field of this 
system is 1.035 and above this, the system exists only in the ferromagnetic state. 

2.4 Dimerized to uniform phase transition 

The dimerized to uniform structural phase transition in the spin-Peierls system 
shown in figure ( |1.3| ) is obtained numerically, for the case of zero magnetic field 
(/i=0) and zero magnetization (m = 0, A'"i7'=50), as discussed below. 

By taking the bond displacements to be of the form 5„ = 1)"", 
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Figure 2.9: Dimerized to uniform phase transition at zero magnetic field. 



which is true for dimerized phase, the value of 5 that minimizes the free energy 
is obtained for each non-zero temperature, by following the same procedure 
mentioned in section (2.2), except that free energy is minimized instead of the 
ground-state energy. The free energy U at any temperature T, using the Fermi- 
Dirac statistics, is given by 



V = -TT^ E log 1 + 2.8 

where E{n) are the single-particle energies at temperature T, and is the 
Boltzmann constant [0. 

The step sizes chosen in varying the temperature and the magnitude 
of the bond displacement while minimizing the free energy are both equal to 
0.0001. Fig. (|2.9|) shows the second order phase transition, the critical param- 
eter being the magnitude of the bond displacement, 6, plotted as a function of 



Chapter 2. Numerical Results 



18 



temperature, T. The critical temperature as obtained from the graph is 0.036. 
It is to be noted that though the system is ID in the Hamiltonian ( |2.6| ) it still 
has a phase transition. This is because our calculation is equivalent to a mean- 
field theory which allows a non-zero displacement 5. In a real CuGeOs system 
the interchain couplings justify this mean- field treatment and led to a finite 
temperature phase transition. 



2.5 Calculation of the phonon frequencies 
2.5.1 Analytical method 

In Section (2.2), single-particle energies were calculated assuming static lattice 
displacements. In this section, we will calculate the phonon frequencies using 
perturbation theory. 



The simplified spin-Peierls Hamiltonian is given by the Eq. (O). 
The i^^ bond displacement, 5j , i, is given by 



-A (wi+i - Ui) 



Now the first order change in the ground-state energy (at zero tem- 
perature) Eq due to the displacement of z*'^ bond proportional to OEq/ dS^j^i. 
This is zero since Eq is calculated by minimizing H in Eq. ( |2.6| ) with respect to 
all the bond displacements. Hence to realize the contribution of bond displace- 
ments to a change in the ground- state energy, we have to find the derivative of 
second order with respect to the bond displacements. 

Defining the elements of the dynamical matrix D to be 



where all the derivatives are evaluated at those values of 5's which are obtained 
by minimizing the Hamiltonian ( p.6| ) 0. Using this definition for D we get the 
lattice Hamiltonian to be 
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H,H = \MY.<^.-^ + \Y.D.H^,^[Ku^+i-ud + 5^.]^^^^^ (2.9) 

where Ui is the displacement of the ion at the lattice point i, M its mass, 6')'^i 
is the equilibrium value of 6^^i which minimizes the Hamiltonian and A is the 
positive constant which is a measure of the coupling that exists between the 
lattice points and spins. 

Now we define a transformation on displacements of lattice points as 
shown below, 



Ui = Ui - Ci (2.10) 

we see that with periodic boundary conditions, and with the choice ci = we 
have 

A(q+i-q) + 55°), =0 

since I]i ^j+i = 0. 

Writing the lattice Hamiltonian ( |2.9| ) in terms of ttj's, thus simplifies 
it to the form 



1 ^ 2 

Hph = +yY1 Dij{ui+i - Ui) {uj+i - Uj) 

Simplifying it further by defining 

A. = A-l.-i + D,^y^ - D^-y^ - A+i.-i (2.11) 
we get the lattice Hamiltonian Hph (|2.9| ) to be 

Hph = ^MJ2ii^' + ^Y.D^,JU,UJ (2.12) 
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Now the frequency, tUi, of the z*'* phonon, (proportional to its energy) 

is given by 




(2.13) 



where jjn is the i eigenvalue of the transformed dynamical matrix, D, and M 
is the mass of each lattice point 1^ . 




10 20 30 40 50 60 70 80 90 100 

i > 

Figure 2.10: Square root of eigenvalues of the transformed dynamical matrix 
for the system in dimerized phase. 

The transformed dynamical matrix for the dimerized phase is eval- 
uated assuming bond displacements of the form 5i = 5(— 1)* and the root of 
the eigenvalues of the corresponding Hamiltonian are plotted in an increasing 
order as shown in Fig. ( p.lOp . This curve gives the phonon dispersion for the 
dimerized phase accurate upto a constant factor given by Eq. ( p.l3|) . While 
evaluating the elements of the transformed dynamical matrix, various symme- 
tries like inversion symmetry about the bonds, translational symmetry of lattice 
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points, etc., are used to generate all the elements from just the first few rows or 
columns. 



2.5.2 Numerical method 



The eigenf unctions of the Hamiltonian, H for the case of dimerized phase (zero 
magnetization) are used to evaluate the transformed dynamical matrix using 
the second order perturbation theory as shown below 



D{m, n) 

< n+ l\H\n >< n\H\n + 1 > + < m + l\H\m >< m|i7|m + 1 > 



E 



m - E{k) 



(2.14) 



where N is the system size, j represents an occupied single-particle state and 
takes the values from 1 to 50, k represents an unoccupied single-particle state 
and ranges from 51 to 100, E[j) and E[k) denote the single-particle energies, 
m and n denote the site numbers with \m > representing a state in which the 
m*'* site is occupied. While most of the matrix elements of D obtained using 
the above equation are found to agree very well with those determined by the 
analytical method described in the previous section, a few off-diagonal elements 
are seen to differ shghtly. 



Chapter 3 
Analytical Methods 



3.1 Uniform phase 

The Hamiltonian for the system in uniform phase is given by 



H = -— [c+c„+i + c:^+^Cn + (-1)^^ c^Ci + c^Cn 

^ n=l 

+ (3.1) 

n=l 



This is obtained by setting each and also the total magnetiza- 
tion, m, to zero in (|2.6| ). The single-particle energy spectrum of this Hamiltonian 



is obtained by taking the discrete Fourier transform for this equation [g, [IG 
This is given by = —Jcos{k), where k is the Fourier momentum and is shown 
in the left half of Fig. ( |3.1D . In the plot, the value of exchange integral J is 
chosen to be unity. Also the lattice constant I of uniform phase is set to unity all 
throughout this report. Imposing periodic boundary conditions on the system, 
k assumes the form k = 27rn/N where N is the system size and n takes values 
from to iV - 1. 

The ground-state energy of the system per site at zero temperature, 
Ek is now obtained by summing the lower half of the energy spectrum and 
dividing it by the system size A^, and is found to be — J/vr. 



22 



Chapter 3. Analytical Methods 



23 



Uniform Phase 



Dimerized Phase 





momentum (k) 



momentum (k) 



Figure 3.1: The single-particle energy spectra for uniform and dimerised phases. 

3.2 Dimerized phase 



In dimerized phase, the bond displacements are of the form = 
(— The total magnetization of the system remains zero. The Hamiltonian 
of the system is shown in (|3.2| ) . 



H 



^ n=l 
^ n=l 



+ (-1 



(3.2) 



The Hamiltonian can be solved using the translation symmetry by 
redefining the size of unit cell in the dimerized phase to be twice that of the 
uniform phase and then applying the continuous Fourier transform as shown be- 
low. In equation ( |3.3| ), A^, B/^ represents the annihilation operators for fermions 
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with momentum k associated with first and second of the pair of lattice points 
present in the unit cell respectively and their Hermitian conjugates represent 
the respective creation operators. 



H 



dk 

W) 



At Bt 



J+S 



J+5 J—& ik 

2 2 ^ 



J—5 —ik 

2 2 ^ 









. Bk _ 



(3.3) 



The single-particle energy spectrum Ck, of this Hamiltonian is given by 



e. = ±i^^^l7l + -^^cosA. 



(3.4) 



where k goes from — vr to vr as shown in right half of Fig. ( p.l|) . In the plot, 
exchange integral J is set to unity. The ground-state energy of the system in 
this dimerized phase at zero temperature is calculated in the following section. 



3.2.1 Ground-state energy 

The total energy of the ground-state is obtained by summing all the negative 
energy states. From last section, it can be seen that the crystal momentum k 
takes values from — tt to vr. Hence, the Ground-state energy per site, Ed is given 
by 



Jl) riT 

Ed = ^ / dky/l + a cos k 

47r V2 J-n 

where a = 1 - 

J2 

6=1 + 1; (3.5) 

This is an elliptical integral of the second kind. Using the series approximation 
for this integral and truncating the series to the quadratic power of [6/ J), 
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^ -^[l-^(i-ai + 2Wog(j))] 

where e(l - ^) = 1 + a,{^) - 26i(^)log(^) 

and ai ^ 0.46301, 6i ^ 0.24527 (3.6) 



Thus, the ground-state energy of the system in dimerized phase is 
lower compared to that of uniform phase for sufficiently smaller values of 6/ J 
and hence dimerization occurs. 



3.3 Modeling the bond displacements for the 
case of magnetization=49 (Np = 99) 

The Hamiltonian for the spin-Peierls system at zero temperature is 
shown in ( |2.6| ). In the continuum limit, the eigenfunctions of this are of the form 
= e*"*^'^^'^'* where the lattice spacing is chosen to be unity and the magnetic 
field, h, is taken to be zero. Substituting this in the above Hamiltonian, we get 



Eo = -(l + (5„+i)(e'('=+") + e-'('=+")) 
= -(l + (5„+i)cos/c 

^ -l+v{x) + — (3.7) 

where Eq is the energy eigenvalue corresponding to ipn and v{x) = —S^_^_i. It 
is to be noted that in the following discussion the lattice label n is replaced by 
the continuum label x whenever convenient. 

Using this form of in the Hamiltonian (|2.6| ) we get 



H = 



(3.8) 
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subject to the constraint / dxv{x) = 0. 

The minimum energy solution of this Hamihonian, ( |3.8|) ,is obtained 
using a trial function ip^x) = 1/ cosh(aa;) where a > 0. Consequently, it reduces 
to 



en —11, , v(x) 

= ( : ) H — — 

cosh(Q;x) cosh(ax) 2 cosh(ax) cosh'^(ax) cosh(ax) 

Hence we choose v{x) = —o?l cosh^(ax). Then matching the coeffi- 
cient of ^/'(a;) = 1/ cosh(aa;) on the two sides gives the equation 



dt 



2 2 J-oo cosh^ t 



(3.10) 



The minimum energy of Cq is attained when a = ^ and it is equal to —1 — 
and the corresponding v{x) = — = — l/4c^ cosh^(^). For c = 1, this value 
is equal to -1.125. From section (2.2.3), numerically this value is found to be 
-1.091. This discrepancy in the value of v(x) arises because the continuum 
approximation is valid only for c >> 1. Also it can be noted that 



dx 



4c2 cosh2(^) 



(3.11) 



So in a finite system with N sites where >> c, we introduce a Lagrangian 
multiplier in Eq. (3^ to enforce the constraint dx J v{x) = 0. We then find 
= — for n » c instead of = O.Then 



eo 



+ 



24c2 2A^c2' 



(3.12) 



Chapter 3. Analytical Methods 



27 



It can be verified that this solution for il){x) is indeed the minimum energy 
solution since the first order change in Hamiltonian ( ^.11 ) is zero. The analytical 
function of obtained here agrees perfectly well with the lower curve in Fig. 
for the values c = 1, iV = 100. 



Chapter 4 



Scope for future work 



4.1 Finite temperature calculations 

All the ground-state energies and eigenfunctions for the magnetiza- 
tions discussed in section (2.2) are calculated at zero temperature. An attempt is 
made to extend these calculations using the same algorithm to non-zero temper- 
atures using Fermi statistics as shown below. The free-energy U at temperature 
T is now given by 

U = -T^ log [1 + exp( )] (4.1) 

where E{n) are the single particle energies at temperature T obtained after min- 
imizing the Hamiltonian ( |2.6| ) and is the Boltzmann constant [0]. However, 
it is noticed that convergence in finding the set of lattice displacements which 
minimize the free-energy could not be achieved. Also, the same problem was 
encountered while calculating the single particle energies and eigenfunctions for 
low-lying excited states at zero temperature. A Solution is yet to be found to 
overcome this problem. This is required to obtain the complete phase diagram 
and compare it with Fig. (|1.3| ). 
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4.2 Calculation of change in eigenfunctions due 
to dynamical phonons 

In section (2.2), a procedure has been detailed to evaluate all the eigenfunc- 
tions assuming the lattice points to be static. The phonon frequencies are then 
calculated in section (2.5) to leading order in lattice displacements by using 
perturbation theory. The corrections to the unperturbed eigenfunctions due to 
these phonons can thereafter be obtained self-consistently. Work to accomplish 
this is in progress. This will give a more accurate picture of the soliton profile in 
the presence of dynamical phonons. For instance, the soliton profiles shown in 
Figs. ( p.3|) and ( p.6|) which do not take phonons into account now get modified. 

4.3 Modeling the lattice spin displacements for 
all possible magnetizations 

The lattice spin displacements for all possible magnetizations can be modelled 
by adopting a procedure similar to the one discussed in section (3.3). A trial 
function for the lattice displacements is used to minimize the Hamiltonian for 
that particular magnetization. The necessary constraints required to be satisfied 
are incorporated in this equation using Lagrangian multipliers. By comparing 
the analytical shapes of the bond displacements with those observed numerically, 
one can find the value of the lattice-spin coupling parameter A. 

4.4 Inclusion of the S^S^^i term in the Hamil- 
tonian 

As detailed in section (2.2.1), for all calculations in this report the generalised 
Heisenberg spins in the full Hamiltonian ( |1.1| ) are replaced by the XY spins. 
This is done to simplify the calculations without qualitatively affecting the phase 
diagram. However, a more realistic picture requires the inclusion of the spin- 
spin interactions in the z direction which on Jordan- Wignerization become four 
fermion interactions 0, Techniques like the density matrix renormailsation 
group and bosonization can be employed to study the Hamiltonian of such a 
system [|]. 
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